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The Gamma distribution has the probability den- 
sity given by a product of an exponential function and a 
power-law function, as in Eq. ([9]). 

The Pareto distribution has the probability density 
P(x) oc l/x 1+a and the cumulative probability C(x) oc 
1 jx a given by a power law. These expressions apply only 
for high enough values of x and do not apply for x — > 0. 

The Lorenz curve was introduced by the American 
economist Max Lorenz to describe income and wealth in- 
equality. It is defined in terms of two coordinates x(r) 
and y(r) given by Eq. (|19[) . The horizontal coordinate 
x(r) is the fraction of the population with income below 
r, and the vertical coordinate y(r) is the fraction of in- 
come this population accounts for. As r changes from 
to oo, x and y change from to 1, parametrically defining 
a curve in the (x, y)-plane. 

The Gini coefficient G was introduced by the Italian 
statistician Corrado Gini as a measure of inequality in a 
society. It is defined as the area between the Lorenz curve 
and the straight diagonal line, divided by the area of the 
triangle beneath the diagonal line. For perfect equality 
(everybody has the same income or wealth) G — 0, and 
for total inequality (one person has all income or wealth, 
and the rest have nothing) G = 1. 

The Fokker-Planck equation is the partial differ- 
ential equation that describes evolution in time t of 
the probability density P(r, t) of a random variable r ex- 
periencing small random changes Ar during short time 
intervals At. It is also known in mathematical litera- 
ture as the Kolmogorov forward equation. The diffusion 
equation is an example of the Fokker-Planck equation. 



Glossary 

Probability density P(x) is defined so that the prob- 
ability of finding a random variable x in the interval from 
x to x + dx is equal to P(x) dx. 

Cumulative probability C(x) is defined as the in- 
tegral C{x) = J™ P(x) dx. It gives the probability that 
the random variable exceeds a given value x. 

The Boltzmann-Gibbs distribution gives the 
probability of finding a physical system in a state with 
the energy e. Its probability density is given by the ex- 
ponential function ([1}. 



I. DEFINITION OF THE SUBJECT 

Econophysics is an interdisciplinary research field ap- 
plying methods of statistical physics to problems in eco- 
nomics and finance. The term "econophysics" was first 
introduced by the prominent theoretical physicist Eugene 
Stanley in 1995 at the conference Dynamics of Com- 
plex Systems, which was held in Calcutta (now known 
as Kolkata) as a satellite meeting to the STATPHYS-19 
conference in China [3, ■ The term appeared in print 
for the first time in the paper by Stanley et al. [3[ in 
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the proceedings of the Calcutta conference. The paper 
presented a manifesto of the new field, arguing that "be- 
havior of large numbers of humans (as measured, e.g., by 
economic indices) might conform to analogs of the scaling 
laws that have proved useful in describing systems com- 
posed of large numbers of inanimate objects" Soon 
the first econophysics conferences were organized: Inter- 
national Workshop on Econophysics, Budapest, 1997 and 
International Workshop on Econophysics and Statistical 
Finance, Palermo, 1998 and the book An Introduc- 
tion to Econophysics [i[ was published. 

The term "econophysics" was introduced by analogy 
with similar terms, such as astrophysics, geophysics, and 
biophysics, which describe applications of physics to dif- 
ferent fields. Particularly important is the parallel with 
biophysics, which studies living creatures, which still 
obey the laws of physics. It should be emphasized that 
econophysics does not literally apply the laws of physics, 
such as Newton's laws or quantum mechanics, to humans, 
but rather uses mathematical methods developed in sta- 
tistical physics to study statistical properties of complex 
economic systems consisting of a large number of hu- 
mans. So, it may be considered as a branch of applied 
theory of probabilities. However, statistical physics is 
distinctly different from mathematical statistics in its fo- 
cus, methods, and results. 

Originating from physics as a quantitative science, 
econophysics emphasizes quantitative analysis of large 
amounts of economic and financial data, which became 
increasingly available with the massive introduction of 
computers and the Internet. Econophysics distances it- 
self from the verbose, narrative, and ideological style of 
political economy and is closer to econometrics in its 
focus. Studying mathematical models of a large number 
of interacting economic agents, econophysics has much 
common ground with the agent-based modeling and 
simulation. Correspondingly, it distances itself from the 
representative-agent approach of traditional economics, 
which, by definition, ignores statistical and heteroge- 
neous aspects of the economy. 

Two major directions in econophysics are applications 
to finance and economics. Applications to finance are de- 
scribed in a separate article, Econophysics of Finan- 
cial Markets, in this encyclopedia. Observational as- 
pects are covered in the article Econophysics, Obser- 
vational. The present article, Econophysics, Statis- 
tical Mechanics Approach to, concentrates primarily 
on statistical distributions of money, wealth, and income 
among interacting economic agents. 

Another direction related to econophysics has been ad- 
vocated by the theoretical physicist Serge Galam since 
early 1980 under the name of sociophysics [a, with the 
first appearance of the term in print in Ref . [y] . It echoes 
the term "physique sociale" proposed in the nineteenth 
century by Auguste Comte, the founder of sociology. Un- 
like econophysics, the term "sociophysics" did not catch 
on when first introduced, but it is coming back with the 
popularity of econophysics and active promotion by some 



physicists @, H, Q- While the principles of both fields 
have a lot in common, econophysics focuses on the nar- 
rower subject of economic behavior of humans, where 
more quantitative data is available, whereas sociophysics 
studies a broader range of social issues. The boundary 
between econophysics and sociophysics is not sharp, and 
the two fields enjoy a good rapport [l(J. A more detailed 
description of historical development in presented in Sec. 

IE 

II. HISTORICAL INTRODUCTION 

Statistical mechanics was developed in the second half 
of the nineteenth century by James Clerk Maxwell, Lud- 
wig Boltzmann, and Josiah Willard Gibbs. These physi- 
cists believed in the existence of atoms and developed 
mathematical methods for describing their statistical 
properties, such as the probability distribution of veloci- 
ties of molecules in a gas (the Maxwell-Boltzmann distri- 
bution) and the general probability distribution of states 
with different energies (the Boltzmann-Gibbs distribu- 
tion). There are interesting connections between the de- 
velopment of statistical physics and statistics of social 
phenomena, which were recently brought up by the sci- 
ence journalist Philip Ball [TJ LI2] • 

Collection and study of "social numbers" , such as the 
rates of death, birth, and marriage, has been growing 
progressively since the seventeenth century [12, Ch. 3]. 
The term "statistics" was introduced in the eighteenth 
century to denote these studies dealing with the civil 
"states", and its practitioners were called "statists". 
Popularization of social statistics in the nineteenth cen- 
tury is particularly accredited to the Belgian astronomer 
Adolphe Quetelet. Before the 1850s, statistics was con- 
sidered an empirical arm of political economy, but then 
it started to transform into a general method of quanti- 
tative analysis suitable for all disciplines. It stimulated 
physicists to develop statistical mechanics in the second 
half of the nineteenth century. 

Rudolf Clausius started development of the kinetic the- 
ory of gases, but it was James Clerk Maxwell who made 
a decisive step of deriving the probability distribution 
of velocities of molecules in a gas. Historical studies 
show 12, Ch. 3] that, in developing statistical mechan- 
ics, Maxwell was strongly influenced and encouraged by 
the widespread popularity of social statistics at the time. 
This approach was further developed by Ludwig Boltz- 
mann, who was very explicit about its origins [12|, p. 69]: 

"The molecules are like individuals, . . . and 
the properties of gases only remain unaltered, 
because the number of these molecules, which 
on the average have a given state, is con- 
stant." 

In his book Populare Schrifen from 1905 [l3j], Boltzmann 
praises Josiah Willard Gibbs for systematic development 
of statistical mechanics. Then, Boltzmann says (cited 
from [3): 
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"This opens a broad perspective, if we do not 
only think of mechanical objects. Let's con- 
sider to apply this method to the statistics of 
living beings, society sociology and so forth." 

(The author is grateful to Michael E. Fisher for bringing 
this quote to his attention.) 

It is worth noting that many now-famous economists 
were originally educated in physics and engineering. Vil- 
fredo Pareto earned a degree in mathematical sciences 
and a doctorate in engineering. Working as a civil engi- 
neer, he collected statistics demonstrating that distribu- 
tions of income and wealth in a society follow a power law 
(l5j |. He later became a professor of economics at Lau- 
sanne, where he replaced Leon Walras, also an engineer 
by education. The influential American economist Irv- 
ing Fisher was a student of Gibbs. However, most of the 
mathematical apparatus transferred to economics from 
physics was that of Newtonian mechanics and classical 
thermodynamics [l6l |. It culminated in the neoclassical 
concept of mechanistic equilibrium where the "forces" 
of supply and demand balance each other. The more 
general concept of statistical equilibrium largely eluded 
mainstream economics. 

With time, both physics and economics became more 
formal and rigid in their specializations, and the social 
origin of statistical physics was forgotten. The situation 
is well summarized by Philip Ball [12, p. 69]: 

"Today physicists regard the application of 
statistical mechanics to social phenomena as 
a new and risky venture. Few, it seems, re- 
call how the process originated the other way 
around, in the days when physical science 
and social science were the twin siblings of a 
mechanistic philosophy and when it was not 
in the least disreputable to invoke the habits 
of people to explain the habits of inanimate 
particles." 

Some physicists and economists attempted to connect 
the two disciplines during the twentieth century. The 
theoretical physicist Ettore Majorana argued in favor of 
applying the laws of statistical physics to social phenom- 
ena in a paper published after his mysterious disappear- 
ance [l7{ ■ The statistical physicist Elliott Montroll co- 
authored the book Introduction to Quantitative Aspects 
of Social Phenomena [l8j . Several economists app lied 
statistical physics to economic problems fl9l 12(1 [2lT[22|. 
An early attempt to bring together the leading theoret- 
ical physicists and economists at the Santa Fe Institute 
was not entirely successful [23]. However, by the late 
1990s, the attempts to apply statistical physics to social 
phenomena finally coalesced into the robust movements 
of econophysics and sociophysics, as described in Sec. |U 

The current standing of econophysics within the 
physics and economics communities is mixed. Although 
an entry on econophysics has a ppe ared in the New Pal- 
grave Dictionary of Economics [241 ] , it is fair to say that 



econophysics is not accepted yet by mainstream eco- 
nomics. Nevertheless, a number of open-minded, nontra- 
ditional economists have joined this movement, and the 
number is growing. Under these circumstances, econo- 
physicists have most of their papers published in physics 
journals. The journal Physica A: Statistical Mechanics 
and its Applications emerged as the leader in econo- 
physics publications and has even attracted submissions 
from some bona fide economists. The mainstream physics 
community is generally sympathetic to econophysics, al- 
though it is not uncommon for econophysics papers to 
be rejected by Physical Review Letters on the grounds 
that "it is not physics". There are regular conference 
in econophysics, such as Applications of Physics in Fi- 
nancial Analysis (sponsored by the European Physical 
Society), Nikkei Econophysics Symposium, and Econo- 
physics Colloquium. Econophysics sessions are included 
in the annual meetings of physical societies and statis- 
tical physics conferences. The overlap with economists 
is the strongest in the field of agent-based simulation. 
Not surprisingly, the conference series WEHIA/ESHIA, 
which deals with heterogeneous interacting agents, regu- 
larly includes sessions on econophysics. 



III. STATISTICAL MECHANICS OF MONEY 
DISTRIBUTION 

When modern econophysics started in the middle of 
1990s, its attention was primarily focused on analysis 
of financial markets. However, three influential papers 
[HI, US H3] ' dealing with the subject of money and wealth 
distributions, were published in year 2000. They started 
a new direction that is closer to economics than finance 
and created an expanding wave of follow-u p pu blications. 
We start reviewing this subject with Ref. [25|], whose re- 
sults are the most closely related to the traditional sta- 
tistical mechanics in physics. 



A. The Boltzmann-Gibbs distribution of energy 

The fundamental law of equilibrium statistical me- 
chanics is the Boltzmann-Gibbs distribution. It states 
that the probability P(e) of finding a physical system or 
sub-system in a state with the energy e is given by the 
exponential function 



P(s) 



ce 



-e/T 



(1) 



where T is the temperature, and c is a normalizing con- 
stant 28]. Here we set the Boltzmann constant ^ to 
unity by choosing the energy units for measuring the 
physical temperature T. Then, the expectation value of 
any physical variable x can be obtained as 



(x) 



(2) 
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where the sum is taken over all states of the system. 
Temperature is equal to the average energy per particle: 
T ~ (e), up to a numerical coefficient of the order of 1. 

Eq. (fTJ) can be derived in different ways [28| . All deriva- 
tions involve the two main ingredients: statistical char- 
acter of the system and conservation of energy s. One 
of the shortest derivations can be summarized as follows. 
Let us divide the system into two (generally unequal) 
parts. Then, the total energy is the sum of the parts: 
e = E\ + £2, whereas the probability is the product of 
probabilities: P(e) — P(ei) Pfa). The only solution of 
these two equations is the exponential function ([T]). 

A more sophisticated derivation, proposed by Boltz- 
mann himself, uses the concept of entropy. Let us con- 
sider N particles with the total energy E. Let us divide 
the energy axis into small intervals (bins) of width Ae 
and count the number of particles Nk having the ener- 
gies from Ek to Ek + Ae. The ratio Nk/N = Pk gives the 
probability for a particle to have the energy Ek- Let us 
now calculate the multiplicity W, which is the number 
of permutations of the particles between different energy 
bins such that the occupation numbers of the bins do 
not change. This quantity is given by the combinatorial 
formula in terms of the factorials 

N\ 

W= — . (3) 

N!\N 2 lNa\ •■■ 

The logarithm of multiplicity is called the entropy S — 
In W . In the limit of large numbers, the entropy per 
particle can be written in the following form using the 
Stirling approximation for the factorials 

k v ' k 

Now we would like to find what distribution of particles 
between different energy states has the highest entropy, 
i.e., the highest multiplicity, provided that the total en- 
ergy of the system, E = J2k^k£k, has a fixed value. 
Solution of this problem can be easily obtained using the 
method of Lagrange multipliers [HI , and the answer gives 
the exponential distribution (fT]). 

The same result can be derived from the ergodic the- 
ory, which says that the many-body system occupies all 
possible states of a given total energy with equal proba- 
bilities. Then it is straightforward to show [29|, |30] that 
the probability distribution of the energy of an individual 
particle is given by Eq. |T]). 



B. Conservation of money 

The derivations outlined in Sec. IIII.Al are very general 
and use only the statistical character of the system and 
the conservation of energy. So, one may expect that the 
exponential Boltzmann-Gibbs distribution ([T]) may apply 
to other statistical systems with a conserved quantity. 



The economy is a big statistical system with millions 
of participating agents, so it is a promising target for ap- 
plications of statistical mechanics. Is there a conserved 
quantity in the economy? The paper [25| argued that 
such a conserved quantity is money m. Indeed, the 
ordinary economic agents can only receive money from 
and give money to other agents. They are not permit- 
ted to "manufacture" money, e.g., to print dollar bills. 
When one agent i pays money Am to another agent j 
for some goods or services, the money balances of the 
agents change as follows 

m.j — > = rrii — Am, 

rrij — > rrij = rrij + Am. (5) 

The total amount of money of the two agents before and 
after transaction remains the same 

nii + rrij = m ■ + m'j , (6) 

i.e., there is a local conservation law for money. The rule 
(JSJ) for the transfer of money is analogous to the transfer 
of energy from one molecule to another in molecular col- 
lisions in a gas, and Eq. @ is analogous to conservation 
of energy in such collisions. 

Addressing some misunderstandings developed in eco- 
nomic literature j^H H3, l33l . [3~i| . we should emphasize 
that, in the model of Ref. [2||, the transfer of money 
from one agent to another happens voluntarily, as a pay- 
ment for goods and services in a market economy. How- 
ever, the model only keeps track of money flow, but does 
not keep track of what kind of goods and service are de- 
livered. One reason for this is that many goods, e.g., 
food and other supplies, and most services, e.g., getting 
a haircut or going to a movie, are not tangible and disap- 
pear after consumption. Because they are not conserved, 
and also because they are measured in different physical 
units, it is not very practical to keep track of them. In 
contrast, money is measured in the same unit (within a 
given country with a single currency) and is conserved 
in transactions, so it is straightforward to keep track of 
money flow. 

Unlike, ordinary economic agents, a central bank or a 
central government can inject money into the economy. 
This process is analogous to an influx of energy into a sys- 
tem from external sources, e.g., the Earth receives energy 
from the Sun. Dealing with these situations, physicists 
start with an idealization of a closed system in thermal 
equilibrium and then generalize to an open system sub- 
ject to an energy flux. As long as the rate of money 
influx from central sources is slow compared with relax- 
ation processes in the economy and does not cause hy- 
perinflation, the system is in quasi-stationary statistical 
equilibrium with slowly changing parameters. This situa- 
tion is analogous to heating a kettle on a gas stove slowly, 
where the kettle has a well-defined, but slowly increasing 
temperature at any moment of time. 

Another potential problem with conservation of money 
is debt. This issue is discussed in more detail in Sec. 
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IIII.DI As a starting point, Ref. first considered sim- 
ple models, where debt is not permitted. This means that 
money balances of agents cannot go below zero: mi > 
for all i. Transaction (0 takes place only when an agent 
has enough money to pay the price: m; > Am, otherwise 
the transaction does not take place. If an agent spends 
all money, the balance drops to zero m t = 0, so the agent 
cannot buy any goods from other agents. However, this 
agent can still produce goods or services, sell them to 
other agents, and receive money for that. In real life, 
cash balance dropping to zero is not at all unusual for 
people who live from paycheck to paycheck. 

The conservation law is the key feature for the success- 
ful functioning of money. If the agents were permitted 
to "manufacture" money, they would be printing money 
and buying all goods for nothing, which would be a dis- 
aster. The physical medium of money is not essential, as 
long as the conservation law is enforced. Money may be 
in the form of paper cash, but today it is more often rep- 
resented by digits in computerized bank accounts. The 
conservation law is the fundamental principle of account- 
ing, whether in the single-entry or the double-entry form. 
More discussion of banks and debt is given in Sec. IIII.DI 

C. The Boltzmann-Gibbs distribution of money 

Having recognized the principle of money conserva- 
tion, Ref. [25[ argued that the stationary distribution of 
money should be given by the exponential Boltzmann- 
Gibbs function analogous to Eq. (jTJ) 

P(m) =ce~ m/T "\ (7) 

Here c is a normalizing constant, and T m is the "money 
temperature" , which is equal to the average amount of 
money per agent: T = (m) = M/N, where M is the total 
money, and N is the number of agents. 

To verify this conjecture, Ref. |25[ performed agent- 
based computer simulations of money transfers between 
agents. Initially all agents were given the same amount 
of money, say, $1000. Then, a pair of agents was 
randomly selected, the amount Am was transferred from 
one agent to another, and the process was repeated many 
times. Time evolution of the probability distribution of 
money P(m) can be seen in computer animation videos at 
the Web pages [H, Hl[ • After a transitory period, money 
distribution converges to the stationary form shown in 
Fig. [TJ As expected, the distribution is very well fitted 
by the exponential function (|7j). 

Several different rules for Am were considered in Ref. 
[25l | . In one model, the transferred amount was fixed 
to a constant Am = $1. Economically, it means that 
all agents were selling their products for the same price 
Am = $1. Computer animation [35| shows that the ini- 
tial distribution of money first broadens to a symmet- 
ric, Gaussian curve, characteristic for a diffusion process. 
Then, the distribution starts to pile up around the m = 
state, which acts as the impenetrable boundary, because 



of the imposed condition m > 0. As a result, Pirn) be- 
comes skewed (asymmetric) and eventually reaches the 
stationary exponential shape, as shown in Fig. [TJ The 
boundary at m = is analogous to the ground state en- 
ergy in statistical physics. Without this boundary con- 
dition, the probability distribution of money would not 
reach a stationary state. Computer animation (35l . [36j 
also shows how the entropy of money distribution, de- 
fined as S/N — — J2k P{ m k) hiP(mfc), grows from the 
initial value S = 0, when all agents have the same money, 
to the maximal value at the statistical equilibrium. 

While the model with Am = 1 is very simple and in- 
structive, it is not very realistic, because all prices are 
taken to be the same. In another model considered in 
Ref. [25j] , Am in each transaction is taken to be a ran- 
dom fraction of the average amount of money per agent, 
i.e., Am = v(M/N), where v is a uniformly distributed 
random number between and 1. The random distri- 
bution of Am is supposed to represent the wide variety 
of prices for different products in the real economy. It 
reflects the fact that agents buy and consume many dif- 
ferent types of products, some of them simple and cheap, 
some sophisticated and expensive. Moreover, different 
agents like to consume these products in different quan- 
tities, so there is variation of paid amounts Am, even 
though the unit price of the same product is constant. 
Computer simulation of this model produces exactly the 
same stationary distribution ([7]), as in the first model. 
Computer animation for this model is also available on 
the Web page [Hj]. 

The final distribution is universal despite different rules 
for Am. To amplify this point further, Ref. [25| also con- 
sidered a toy model, where Am was taken to be a ran- 
dom fraction of the average amount of money of the two 
agents: Am = v[rrii + rrij)/2. This model produced the 
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FIG. 1 Histogram and points: Stationary probability distri- 
bution of money P(m) obtained in agent-based computer sim- 
ulations. Solid curves: Fits to the Boltzmann-Gibbs law ©. 
Vertical lines: The initial distribution of money. (Reproduced 
from Ref. [U) 
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same stationary distribution (JT]) as the two other models. 

The pairwise models of money transfer are attractive 
in their simplicity, but they represent a rather primitive 
market. Modern economy is dominated by big firms, 
which consist of many agents, so Ref. [25| also studied a 
model with firms. One agent at a time is appointed to be- 
come a "firm" . The firm borrows capital K from another 
agent and returns it with interest hK , hires L agents and 
pays them wages W, manufactures Q items of a product, 
sells them to Q agents at price R, and receives profit 
F = RQ — LW — hK. All of these agents are randomly 
selected. The parameters of the model are optimized fol- 
lowing a procedure from economics textbooks (37 ]. The 
aggregate demand-supply curve for the product is taken 
in the form R(Q) = vfQp , where Q is the quantity con- 
sumers would buy at price R, and r\ and v are some 
parameters. The production function of the firm has 
the traditional Cobb-Douglas form: Q(L, K) = L X K 1 ^ X 1 
where x is a parameter. Then the profit of the firm F is 
maximized with respect to K and L. The net result of 
the firm activity is a many-body transfer of money, which 
still satisfies the conservation law. Computer simulation 
of this model generates the same exponential distribution 
0, independently of the model parameters. The reasons 
for the universality of the Boltzmann-Gibbs distribution 
and its limitations are discussed from a different perspec- 
tive in Sec. UlLFl 

Well after the paper [25[ appeared, Italian econophysi- 
cists [38| found that similar ideas had been published 
earlier in obscure journals in Italian by Eleonora Ben- 
nati [H, H3|. They proposed calling these models the 
Bennati-Dragulescu-Yakovenko (BDY) game 41]. The 
Boltzmann distribution was independently applied to so- 
cial sciences by Jiirgen Mimkes using the Lagrange prin- 
ciple of maximization with constraints [42l [43| . The ex- 
ponential distribution of money was also found in Ref. 
(44j using a Markov chain approach to strategic market 
games. A long time ago, Benoit Mandelbrot observed 
M P 83]: 



"There is a great temptation to consider the 
exchanges of money which occur in economic 
interaction as analogous to the exchanges of 
energy which occur in physical shocks be- 
tween gas molecules." 



He realized that this process should result in the expo- 
nential distribution, by analogy with the barometric dis- 
tribution of density in the atmosphere. However, he dis- 
carded this idea, because it does not produce the Pareto 
power law, and proceeded to study the stable Levy distri- 
butions. Ironically, the actual economic data, discussed 
in Sees. IIV.CI and IV. Al do show the exponential distri- 
bution for the majority of the population. Moreover, the 
data have finite variance, so the stable Levy distributions 
are not applicable because of their infinite variance. 



D. Models with debt 

Now let us discuss how the results change when debt is 
permitted. Debt may be considered as negative money. 
When an agent borrows money from a bank (considered 
here as a big reservoir of money), the cash balance of 
the agent (positive money) increases, but the agent also 
acquires a debt obligation (negative money), so the total 
balance (net worth) of the agent remains the same, and 
the conservation law of total money (positive and neg- 
ative) is satisfied. After spending some cash, the agent 
still has the debt obligation, so the money balance of 
the agent becomes negative. Any stable economic sys- 
tem must have a mechanism preventing unlimited bor- 
rowing and unlimited debt. Otherwise, agents can buy 
any products without producing anything in exchange by 
simply going into unlimited debt. The exact mechanisms 
of limiting debt in the real economy are complicated and 
obscured. Ref. 25] considered a simple model where the 
maximal debt of any agent is limited by a certain amount 
nid- This means that the boundary condition m; > is 
now replaced by the condition m, > — md for all agents i. 
Setting interest rates on borrowed money to be zero for 
simplicity, Ref. [25( ] performed computer simulations of 
the models described in Sec. IIII.CI with the new bound- 
ary condition. The results are shown in Fig. O Not 
surprisingly, the stationary money distribution again has 
the exponential shape, but now with the new boundary 
condition at m = —md and the higher money tempera- 
ture Td = rnd + M/N. By allowing agents to go into debt 
up to m,d, we effectively increase the amount of money 
available to each agent by md- So, the money tempera- 
ture, which is equal to the average amount of effectively 
available money per agent, increases. A model with non- 
zero interest rates was also studied in Ref. [25j ]. 

We see that debt does not violate the conservation law 
of money, but rather modifies boundary conditions for 
P(m). When economics textbooks describe how "banks 



create money" or "debt creates money" [37J, they count 
only positive money (cash) as money, but do not count 
liabilities (debt obligations) as negative money. With 
such a definition, money is not conserved. However, if 
we include debt obligations in the definition of money, 
then the conservation law is restored. This approach is in 
agreement with the principles of double-entry accounting, 
which records both assets and debts. Debt obligations 
are as real as positive cash, as many borrowers painfully 
realized in their experience. A more detailed study of 
positive and negative money and book-keeping from the 
point of view of econophysics is presented in a series of 
papers by the physicist Dieter Braun [46|, [47], EH . 

One way of limiting the total debt in the economy is 
the so-called required reserve ratio r (37] | . Every bank is 
required by law to set aside a fraction r of the money 
deposited into the bank, and this reserved money cannot 
be loaned further. If the initial amount of money in the 
system (the money base) is Mo, then with loans and bor- 
rowing the total amount of positive money available to 
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the agents increases to M = Mn/r, where the factor 1/r 
is called the money multiplier (37(. This is how "banks 
create money" . Where does this extra money come from? 
It comes from the increase of the total debt in the sys- 
tem. The maximal total debt is equal to D = Mq/t — Mq 
and is limited by the factor r. When the debt is maxi- 
mal, the total amounts of positive, Mg/r, and negative, 
Mq(1 — r)/r, money circulate between the agents in the 
system, so there arc effectively two conservation laws for 
each of them (49|. Thus, we expect to see the exponential 
distributions of positive and negative money character- 
ized by two different temperatures: T + = Mo/rN and 
T_ = Mq(1 — r)/rN. This is exactly what was found 
in computer simulations in Ref. (49l |. shown in Fig. [3] 
Similar two-sided distributions were also found in Ref. 

E. Proportional money transfers and saving propensity 

In the models of money transfer considered thus far, 
the transferred amount Am is typically independent of 
the money balances of agents. A different model was in- 
troduced in physics literature earlier [50l | under the name 
multiplicative asset exchange model. This model also sat- 
isfies the conservation law, but the transferred amount of 
money is a fixed fraction 7 of the payer's money in Eq. 
©: ' 

Am = 7m,. (8) 

The stationary distribution of money in this model, 
shown in Fig. 2] with an exponential function, is close, 
but not exactly equal, to the Gamma distribution: 

P(m)=cm' 3 e- m/T . (9) 



Eq. ([9]) differs from Eq. ([7]) by the power-law prefactor 
mr . From the Boltzmann kinetic equation (discussed in 
more detail in Sec. IIII.F[) . Ref. [50] derived a formula 
relating the parameters 7 and (3 in Eqs. © and ©: 

/3 = -l-ln2/ln(l-7). (10) 

When payers spend a relatively small fraction of their 
money 7 < 1/2, Eq. (|10p gives (3 > 0, so the low-money 
population is reduced and P(m — > 0) = 0, as shown in 
Fig.rj 

Later, the economist Thomas Lux brought to the 
attention of physicists [32| that essentially the same 
model, called the inequality process, had been introduced 
and studied much earlier by the sociologist John Angle 
[5ll IH2I [H, Hi], HH| , see also the review for additional 
references. While Ref. [5(| did not give much justification 
for the proportionality law ([5]) , Angle [5l| connected this 
rule with the surplus theory of social stratification [57| . 
which argues that inequality in human society develops 
when people can produce more than necessary for mini- 
mal subsistence. This additional wealth (surplus) can be 
transferred from original producers to other people, thus 
generating inequality. In the first paper by Angle [5l| . 
the parameter 7 was randomly distributed, and another 
parameter S gave a higher probability of winning to the 
agent with a higher money balance in Eq. ([5]). However, 
in the following papers, he simplified the model to a fixed 
7 (denoted as to by Angle) and equal probabilities of win- 
ning for higher- and lower-balance agents, which makes 
it completely equivalent to the model of Ref. [50] • Angle 
also considered a model [H, [56| where groups of agents 
have different values of 7, simulating the effect of edu- 
cation and other "human capital" . All of these models 
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FIG. 2 Histograms: Stationary distributions of money with 
and without debt. The debt is limited to = 800. Solid 
curves: Fits to the Boltzmann-Gibbs laws with the "temper- 
atures" T = 1800 and T = 1000. (Reproduced from Ref. 

M) 
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FIG. 3 The stationary distribution of money for the required 
reserve ratio r = 0.8. The distribution is exponential for 
positive and negative money with different "temperatures" 
T+ and T_, as illustrated by the inset on log-linear scale. 
(Reproduced from Ref. [Is|) 
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generate a Gamma-like distribution, well approximated 
by Eq. ©. 

Another model with an element of proportionality was 
proposed in Ref. [2(| . (This paper originally appeared as 
a follow-up preprint cond-mat/0004256 to the preprint 

this 



cond-mat/0001432| of Ref. |25|J.) In - this model, the 
agents set aside (save) some fraction of their money Amj, 
whereas the rest of their money balance (1 — A)mj be- 
comes available for random exchanges. Thus, the rule of 
exchange ([5]) becomes 



Am, + £(1 - A)(m, + m 3 ), 
: Xmj + (1 - 0(1 - A) (m^ -+ 



(11) 



Here the coefficient A is called the saving propensity, and 
the random variable £ is uniformly distributed between 
and 1. It was pointed out in Ref. [5a| that, by the change 
of notation A — > (1 — 7), Eq. l(TTj) can be transformed to 
the same form as Eq. ((SJ, if the random variable £ takes 
only discrete values and 1. Computer simulations [26j 
of the model (fTTj) found a stationary distribution close 
to the Gamma distribution ©. It was shown that the 
parameter (3 is related to the saving propensity A by the 
formula (3 = 3A/(1 - A) H M, SI For A ^ 0, 

agents always keep some money, so their balances never 
go to zero and P(m — > 0) = 0, whereas for A = the 
distribution becomes exponential. 

In the subsequent papers by the Kolkata school [l[ and 
related papers, the case of random saving propensity was 
studied. In these models, the agents are assigned random 
parameters A drawn from a uniform distribution between 
and 1 [6l|. It was found that this model produces a 
power-law tail P(m) oc 1/m 2 at high m. The reasons 
for stability of this law were understood using the Boltz- 
mann kinetic equation [rjfl . l62l , |63| , but most elegantly in 
the mean- field theory [64|. The fat tail originates from 
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FIG. 4 Histogram: Stationary probability distribution of 
money in the multiplicative random exchange model (JSJ) for 
7 = 1/3. Solid curve: The exponential Boltzmann-Gibbs law. 
(Reproduced from Ref. [25j ~) 



the agents whose saving propensity is close to 1, who 
hoard money and do not give it back [381 [65| . An inter- 
esting matrix formulation of the problem was presented 
in Ref. [6||. Ref. [67j studied the relaxation rate in the 
money transfer models. Ref. [25| studied a model with 
taxation, which also has an element of proportionality. 
The Gamma distribution was also studied for conserva- 
tive models within a simple Boltzmann approach in Ref. 
[68j and using much more complicated rules of exchange 
in Ref. M,M- 



F. Additive versus multiplicative models 

The stationary distribution of money ([9]) for the mod- 
els of Sec. IIII. El is different from the simple exponential 
formula ((7|) found for the models of Sec. IIII. Cl The origin 
of this difference can be understood from the Boltzmann 
kinetic equation [28|, [7l[. This equation describes time 
evolution of the distribution function P(m) due to pair- 
wise interactions: 



dP{m) 
~dt 



{-/[m,m']— [m-A,m'+A]-P(w)P(m / ) (12) 



+f[ m -A,m'+A]^[m,m']P(m - A)P(m' + A)} dm' dA. 

Here f[ m ,m']-+{m-A,m'+A] is the probability of transfer- 
ring money A from an agent with money m to an agent 
with money m! per unit time. This probability, multi- 
plied by the occupation numbers P(m) and P(m'), gives 
the rate of transitions from the state [m, m'\ to the state 
[m - A, rri + A). The first term in Eq. ((T2J) gives the de- 
population rate of the state m. The second term in Eq. 
(fT2|) describes the reversed process, where the occupation 
number P(m) increases. When the two terms are equal, 
the direct and reversed transitions cancel each other sta- 
tistically, and the probability distribution is stationary: 
dP(m)/dt = 0. This is the principle of detailed balance. 

In physics, the fundamental microscopic equations of 
motion of particles obey time-reversal symmetry. This 
means that the probabilities of the direct and reversed 
processes are exactly equal: 



/[ 



7n,m']— *[m — A,m'+A] 



m— A,m'+A] — »[m,m'] • 



(13) 



When Eq. (TIB"]) is satisfied, the detailed balance condi- 
tion for Eq. (TT2"]) reduces to the equation P(m)P(m') — 
P(m — A)P(m' + A), because the factors / cancels out. 
The only solution of this equation is the exponential func- 
tion P(m) = cexp(— m/T m ), so the Boltzmann-Gibbs 
distribution is the stationary solution of the Boltzmann 
kinetic equation (TT2")) . Notice that the transition prob- 
abilities (fl3|) are determined by the dynamical rules of 
the model, but the equilibrium Boltzmann-Gibbs distri- 
bution does not depend on the dynamical rules at all. 
This is the origin of the universality of the Boltzmann- 
Gibbs distribution. It shows that it may be possible to 
find out the stationary distribution without knowing de- 
tails of the dynamical rules (which are rarely known very 
well), as long as the symmetry condition (fl~3"|) is satisfied. 
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The models considered in Sec. IIII.CI have the time- 
reversal symmetry. The model with the fixed money 
transfer A has equal probabilities (fT"3|) of transferring 
money from an agent with balance m to an agent with 
balance m' and vice versa. This is also true when A is 
random, as long as the probability distribution of A is 
independent of m and ml . Thus, the stationary distribu- 
tion P(m) is always exponential in these models. 

However, there is no fundamental reason to expect 
time-reversal symmetry in economics, so Eq. (|13p may 
be not valid. In this case, the system may have a non- 
exponential stationary distribution or no stationary dis- 
tribution at all. In model ([8|) , the time- reversal symmetry 
is broken. Indeed, when an agent i gives a fixed fraction 
7 of his money m, to an agent with balance mj, their 
balances become (1 — 7)771, and rrij +^rrii. If we try to re- 
verse this process and appoint the agent j to be the payer 
and to give the fraction 7 of her money, 7(771.7 + 7^), 
to the agent i, the system does not return to the origi- 
nal configuration [m,, rrij]. As emphasized by Angle [56| . 
the payer pays a deterministic fraction of his money, but 
the receiver receives a random amount from a random 
agent, so their roles are not interchangeable. Because the 
proportional rule typically violates the time-reversal sym- 
metry, the stationary distribution P(m) in multiplicative 
models is typically not exactly exponential. 1 Making the 
transfer dependent on the money balance of the payer 
effectively introduces a Maxwell's demon into the model. 
That is why the stationary distribution is not exponen- 
tial, and, thus, does not maximize entropy (|4]). Another 
view on the time-reversal symmetry in economic dynam- 
ics is presented in Ref. [72j]. 

These examples show that the Boltzmann-Gibbs dis- 
tribution does not hold for any conservative model. How- 
ever, it is universal in a limited sense. For a broad class of 
models that have time-reversal symmetry, the stationary 
distribution is exponential and does not depend on the 
details of the model. Conversely, when the time-reversal 
symmetry is broken, the distribution may depend on the 
details of the model. The difference between these two 
classes of models may be rather subtle. Deviations from 
the Boltzmann-Gibbs law may occur only if the transition 
rates / in Eq. (fl3|) explicitly depend on the agent's money 
m or m! in an asymmetric manner. Ref. [25| performed a 
computer simulation where the direction of payment was 
randomly fixed in advance for every pair of agents 
In this case, money flows along directed links between 
the agents: i —> j —> k, and the time-reversal symme- 
try is strongly violated. This model is closer to the real 
economy, where one typically receives money from an em- 
ployer and pays it to a grocery store. Nevertheless, the 
Boltzmann-Gibbs distribution was found in this model, 



However, when Am is a fraction of the total money mi+rrij of the 
two agents, the model is time-reversible and has an exponential 
distribution, as discussed in Sec. IIII.CI 



because the transition rates / do not explicitly depend 
on m and m' and do not violate Eq. (|13[) . 

In the absence of detailed knowledge of real micro- 
scopic dynamics of economic exchanges, the semiuniver- 
sal Boltzmann-Gibbs distribution ([7]) is a natural starting 
point. Moreover, the assumption of Ref. [25| that agents 
pay the same prices Am for the same products, indepen- 
dent of their money balances m, seems very appropriate 
for the modern anonymous economy, especially for pur- 
chases over the Internet. There is no particular empirical 
evidence for the proportional rules (jHJl or (fTTj) . However, 
the difference between the additive ([7]) and multiplicative 
© distributions may be not so crucial after all. From the 
mathematical point of view, the difference is in the im- 
plementation of the boundary condition at m = 0. In the 
additive models of Sec. IIII.CI there is a sharp cut-off of 
P(m) ^ at 777 = 0. In the multiplicative models of Sec. 
IIII. El the balance of an agent never reaches m = 0, so 
P{m) vanishes at m — ► in a power-law manner. At the 
same time, P{m) decreases exponentially for large m for 
both models. 

By further modifying the rules of money transfer and 
introducing more parameters in the models, one can ob- 
tain even more complicated distributions [73|. However, 
one can argue that parsimony is the virtue of a good 
mathematical model, not the abundance of additional 
assumptions and parameters, whose correspondence to 
reality is hard to verify. 



IV. STATISTICAL MECHANICS OF WEALTH 
DISTRIBUTION 

In the econophysics literature on exchange models, the 
terms "money" and "wealth" are often used interchange- 
ably; however, economists emphasize the difference be- 
tween these two concepts. In this section, we review the 
models of wealth distribution, as opposed to money dis- 
tribution. 



A. Models with a conserved commodity 

What is the difference between money and wealth? On 
can argue [25| that wealth Wi is equal to money m, plus 
the other property that an agent i has. The latter may 
include durable material property, such as houses and 
cars, and financial instruments, such as stocks, bonds, 
and options. Money (paper cash, bank accounts) is gen- 
erally liquid and countable. However, the other property 
is not immediately liquid and has to be sold first (con- 
verted into money) to be used for other purchases. In 
order to estimate the monetary value of property, one 
needs to know the price p. In the simplest model, let us 
consider just one type of property, say, stocks s. Then 
the wealth of an agent i is given by the formula 

Wi = m l +psi. (14) 
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It is assumed that the price p is common for all agents 
and is established by some kind of market process, such 
as an auction, and may change in time. 

It is reasonable to start with a model where both the 
total money M =T- m,: and the total stock S = J2i s « 
are conserved [z3, l75ll7t| . The agents pay money to buy 
stock and sell stock to get money, and so on. Although 
M and S are conserved, the total wealth W — J2i w i 
is g enerally not conserved, because of price fluctuation 
[75| | in Eq. (TH| . This is an important difference from the 
money transfers models of Sec. IIII1 Here the wealth Wi 
of an agent i, not participating in any transactions, may 
change when transactions between other agents establish 
a new price p. Moreover, the wealth Wi of an agent i does 
not change after a transaction with an agent j. Indeed, 
in exchange for paying money Am, agent i receives the 
stock As = Am/p, so her total wealth (TTJ| remains the 
same. In principle, the agent can instantaneously sell the 
stock back at the same price and recover the money paid. 
If the price p never changes, then the wealth Wi of each 
agent remains constant, despite transfers of money and 
stock between agents. 

We see that redistribution of wealth in this model is 
directly related to price fluctuations. One mathematical 
model of this process was studied in Ref. [77j]. In this 
model, the agents randomly change preferences for the 
fraction of their wealth invested in stocks. As a result, 
some agents offer stock for sale and some want to buy 
it. The price p is determined from the market-clearing 
auction matching supply and demand. Ref. [77j demon- 
strated in computer simulations and proved analytically 
using the theory of Markov processes that the stationary 
distribution P(w) of wealth w in this model is given by 
the Gamma distribution, as in Eq. ^ . Various modifica- 
tions of this model [321 ] , such as introducing monopolistic 
coalitions, do not change this result significantly, which 
shows the robustness of the Gamma distribution. For 
models with a conserved commodity, Ref. [75l ] found the 
Gamma distribution for a fixed saving propensity and a 
power law tail for a distributed saving propensity. 

Another model with conserved money and stock was 
studied in Ref. [78j for an artificial stock market, where 
traders follow different investment strategies: random, 
momentum, contrarian, and fundamentalist. Wealth dis- 
tribution in the model with random traders was found 
have a power-law tail P(w) ~ 1/w 2 for large w. How- 
ever, unlike in most other simulation, where all agents 
initially have equal balances, here the initial money and 
stock balances of the agents were randomly populated 
according to a power law with the same exponent. This 
raises the question whether the observed power-law dis- 
tribution of wealth is an artifact of the initial conditions, 
because equilibrization of the upper tail may take a very 
long simulation time. 



B. Models with stochastic growth of wealth 

Although the total wealth W is not exactly conserved 
in the models considered in Sec. IIV.AI nevertheless W 
remains constant on average, because the total money M 
and stock S are conserved. A different model for wealth 
distribution was proposed in Ref. [27|- In this model, 
time evolution of the wealth w; of an agent i is given by 
the stochastic differential equation 



dwi 



Vi(t) m 



Ji i w i 



(15) 



where r\i (t) is a Gaussian random variable with the mean 
(rj) and the variance 2a 2 . This variable represents growth 
or loss of wealth of an agent due to investment in stock 
market. The last two terms describe transfer of wealth 
between different agents, which is taken to be propor- 
tional to the wealth of the payers with the coefficients Jy . 
So, the model (|15p is multiplicative and invariant under 
the scale transformation Wi — * Zwi. For simplicity, the 
exchange fractions are taken to be the same for all agents: 
Jij = J/N for all i ^ j, where N is the total number of 
agents. In this case, the last two terms in Eq. (fT5"|) can 
be written as J((w) — Wi), where (w) = J2i W i/N is the 
average wealth per agent. This case represents a "mean- 
field" model, where all agents feel the same environment. 
It can be easily shown that the average wealth increases 
in time as (w) t — (w}oe"' n ' + ' r Then, it makes more 
sense to consider the relative wealth ibi — Wi/{w)t- Eq. 
(fT5|) for this variable becomes 

^ = W) - (v) t 2 ) Wi + J(l - Wi). (16) 

The probability distribution P(w, t) for the stochastic 
differential equation (I16| is governed by the Fokker- 
Planck equation 

dP _ d[J(w - 1) + cj 2 w]P | j2 d f_ d{wP) \ 
dt dw dw V dw ) 



The stationary solution (dP/dt 
given by the following formula 



0) of this equation is 



P(w) 



, — J/l7 2 W 



(18) 



The distribution (|18p is quite different from the 
Boltzmann-Gibbs (O and Gamma |9]) distributions. Eq. 
(fT5)l has a power-law tail at large w and a sharp cutoff at 
small w. Eq. (|T5|) is a version of the generalized Lotka- 
Volterra model, and the stationary distribution (fTS)) was 
also obtained in Ref. .79, [80j. The model was generalized 
to include negative wealth in Ref. |8l| . 

Ref. [ItJ used the mean-field approach. A similar result 
was found for a model with pairwise interaction between 
agents in Ref. 82]. In this model, wealth is transferred 
between the agents using the proportional rule l[5]). In 
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addition, the wealth of the agents increases by the factor 
1 + £ in each transaction. This factor is supposed to re- 
flect creation of wealth in economic interactions. Because 
the total wealth in the system increases, it makes sense 
to consider the distribution of relative wealth P(w). In 
the limit of continuous trading, Ref. [82| found the same 
stationary distribution (|18p . This result was reproduced 
using a mathematically more involved treatment of this 
model in Ref. [83| . Numerical simulations of the models 
with stochastic noise rj in Ref. [69j, |70j also found a power 
law tail for large w. 

Let us contrast the models discussed in Sees. IIV. Al and 
IIV.B1 In the former case, where money and commodity 
are conserved, and wealth does not grow, the distribu- 
tion of wealth is given by the Gamma distribution with 
the exponential tail for large w. In the latter models, 
wealth grows in time exponentially, and the distribution 
of relative wealth has a power law tail for large w. These 
results suggest that the presence of a power-law tail is 
a nonequilibrium effect that requires constant growth or 
inflation of the economy, but disappears for a closed sys- 
tem with conservation laws. 

Reviews of the discussed models were also given in 
Refs. [13, IH|. Because of lack of space, we omit dis- 
cussion of models with wealth condensation [U H3, HI, 
H3j HI] , where few agents accumulate a finite fraction of 
total wealth, and studies of wealth distribution on net- 
works [1^, Ho, [Sli [52] . In this section, we discussed the 
models with long-range interaction, where any agent can 
exchange money and wealth with any other agent. A 
local model, where agents trade only with the nearest 
neighbors, was studied in Ref. (93j . 



C. Empirical data on money and wealth distributions 

It would be very interesting to compare theoretical re- 
sults for money and wealth distributions in various mod- 
els with empirical data. Unfortunately, such empirical 
data are difficult to find. Unlike income, which is dis- 
cussed in Sec. |V1 wealth is not routinely reported by the 
majority of individuals to the government. However, in 
many countries, when a person dies, all assets must be 
reported for the purpose of inheritance tax. So, in prin- 
ciple, there exist good statistics of wealth distribution 
among dead people, which, of course, is different from the 
wealth distribution among living people. Using an ad- 
justment procedure based on the age, gender, and other 
characteristics of the deceased, the UK tax agency, the 
Inland Revenue, reconstructed the wealth distribution of 
the whole population of the UK (94[ . Fig. [5] shows the 
UK data for 1996 reproduced from Ref. [95|]. The figure 
shows the cumulative probability C(w) = P(w') dw' 
as a function of the personal net wealth w, which is com- 
posed of assets (cash, stocks, property, household goods, 
etc.) and liabilities (mortgages and other debts). Because 
statistical data are usually reported at non-uniform in- 
tervals of w, it is more practical to plot the cumulative 



probability distribution C(w) rather than its derivative, 
the probability density P(w). Fortunately, when P(w) 
is an exponential or a power-law function, then C(w) is 
also an exponential or a power-law function. 

The main panel in Fig. [5] shows a plot of C (w) on a 
log-log scale, where a straight line represents a power-law 
dependence. The figure shows that the distribution fol- 
lows a power law C(w) ex l/w a with exponent a = 1.9 for 
the wealth greater than about 100 k£. The inset in Fig. 
[5] shows the data on log- linear scale, where a straight line 
represents an exponential dependence. We observe that, 
below 100 k£, the data is well fitted by the exponen- 
tial distribution C(w) oc exp(— w/T w ) with the effective 
"wealth temperature" T w = 60 k£ (which corresponds 
to the median wealth of 41 k£). So, the distribution 
of wealth is characterized by the Pareto power law in 
the upper tail of the distribution and the exponential 
Boltzmann-Gibbs law in the lower part of the distribu- 
tion for the great majority (about 90%) of the population. 
Similar results are found for the distribution of income, 
as discussed in Sec. |Vj One may speculate that wealth 
distribution in the lower part is dominated by distribu- 
tion of money, because the corresponding people do not 
have other significant assets, so the results of Sec. IIIII give 
the Boltzmann-Gibbs law. On the other hand, the upper 
tail of wealth distribution is dominated by investment 
assess, where the results of Sec. IIV.BI give the Pareto 
law. The power law was studied by many researchers for 
the upper-tail data, such as the Forbes list of 400 richest 
people [sl, [93 , but much less attention was paid to the 
lower part of the wealth distribution. Curiously, Ref. [98[ 
found that the wealth distribution in the ancient Egyp- 
tian society was consistent with Eq. (fT5|) . 

For direct comparison with the results of Sec. IIIII it 
would be very interesting to find data on the distribu- 



United Kingdom, IR data for 1996 
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FIG. 5 Cumulative probability distribution of net wealth in 
the UK shown on log-log (main panel) and log-linear (inset) 
scales. Points represent the data from the Inland Revenue, 
and solid lines are fits to the exponential (Boltzmann-Gibbs) 
and power (Pareto) laws. (Reproduced from Ref. |95|P 
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tion of money, as opposed to the distribution of wealth. 
Making a reasonable assumption that most people keep 
most of their money in banks, one can approximate the 
distribution of money by the distribution of balances on 
bank accounts. (Balances on all types of bank accounts, 
such as checking, saving, and money manager, associated 
with the same person should be added up.) Despite im- 
perfections (people may have accounts in different banks 
or not keep all their money in banks), the distribution 
of balances on bank accounts would give valuable infor- 
mation about the distribution of money. The data for 
a big enough bank would be representative of the distri- 
bution in the whole economy. Unfortunately, it has not 
been possible to obtain such data thus far, even though 
it would be completely anonymous and not compromise 
privacy of bank clients. 

Measuring the probability distribution of money would 
be very useful, because it determines how much people 
can, in principle, spend on purchases without going into 
debt. This is different from the distribution of wealth, 
where the property component, such as house, car, or re- 
tirement investment, is effectively locked up and, in most 
cases, is not easily available for consumer spending. So, 
although wealth distribution may reflect the distribution 
of economic power, the distribution of money is more rel- 
evant for consumption. Money distribution can be useful 
for determining prices that maximize revenue of a man- 
ufacturer [25(. If a price p is set too high, few people 
can afford it, and, if a price is too low, the sales revenue 
is small, so the optimal price must be in between. The 
fraction of population who can afford to pay the price p 
is given by the cumulative probability C(p), so the total 
sales revenue is proportional to pC (p) . For the exponen- 
tial distribution C(p) = exp(— p/T m ), the maximal rev- 
enue is achieved at p = T m , i.e., the optimal pric e is equal 
to the average amount of money per person [251 ]. Indeed, 
the prices of mass-market consumer products, such as 
computers, electronics, and appliances, remain stable for 
many years at a level determined by their affordability 
to the population, whereas technical parameters of these 
products at the same price level improve dramatically 
owing to technological progress. 



V. DATA AND MODELS FOR INCOME DISTRIBUTION 

In contrast to money and wealth distributions, a lot 
more empirical data are available for the distribution of 
income r from tax agencies and population surveys. In 
this section, we first present empirical data on income 
distribution and then discuss theoretical models. 



A. Empirical data on income distribution 

Empirical studies of income distrib ution hav e a long 
history in the economic literature [9^, llOCt 1 1 ll | . Follow- 
ing the work by Pareto [l5| , much attention was focused 



on the power-law upper tail of the income distribution 
and less on the lower part. In contrast to more compli- 
cated functions discussed in literature, Ref. |l02j | intro- 
duced a new idea by demonstrating that the lower part of 
income distribution can be well fitted with a simple ex- 
ponential function P(r) = cexp(— r/T r ) characterized by 
just one parameter, the "income temperature" T r . Then 
it was recognized that the whole income distribution can 
be fitted by an exponential function in the lowe r part and 
a power-law function in the upper part [95l . ll03l |. as shown 
in Fig. [SI The straight line on the log-linear scale in the 
inset of Fig. [5] demonstrates the exponential Boltzmann- 
Gibbs law, and the straight line on the log-log scale in 
the main panel illustrates the Pareto power law. The fact 
that income distribution consists of two distinct parts 
reveals th e two-class structure of the American society 
[l04 Il05| . Coexistence of the exponential and power- 
law distributions is also known in plasma physics and 
astrophysics, where t hey are called the "thermal" and 
"superthermal" parts [ToS [ToJ Gill ■ The boundary be- 
tween the lower and upper classes can be defined as the 
intersection point of the exponential and power-law fits 
in Fig.0 For 1997, the annual income separating the two 
classes was about 120 k$. About 3% of the population 
belonged to the upper class, and 97% belonged to the 
lower cl ass. 

Ref. [l05| studied time evolution of income distribu- 
tion in the USA during 1983-2001 on the basis of the data 
from the Internal Revenue Service (IRS), the government 
tax agency. The structure of the income distribution was 
found to be qualitatively the same for all years, as shown 
in Fig. [7] The average income in nominal dollars approx- 
imately doubled during this time interval. So, the hori- 
zontal axis in Fig. [7] shows the normalized income r/T r , 
where the "income temperature" T r was obtained by fit- 
ting of the exponential part of the distribution for each 



United States, IRS data for 1 997 




Adjusted Gross Income, k$ 

FIG. 6 Cumulative probability distribution of tax returns for 
USA in 1997 shown on log-log (main panel) and log-linear 
(inset) scales. Points represent the Internal Revenue Service 
data, and solid lines are fits to the expo nential and power-law 
functions. (Reproduced from Ref. |l03f |) 
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year. The values of T r are shown in Fig. [Jj The plots for 
the 1980s and 1990s are shifted vertically for clarity. We 
observe that the data points in the lower-income part of 
the distribution collapse on the same exponential curve 
for all years. This demonstrates that the shape of the in- 
come distribution for the lower class is extremely stable 
and does not change in time, despite gradual increase of 
the average income in nominal dollars. This observation 
suggests that the lower-class distribution is in statistical, 
"thermal" equilibrium. 

On the other hand, Fig. [JJ shows that the income dis- 
tribution in the upper class d oes n ot rescale and signifi- 
cantly changes in time. Ref. |l05l ] found that the expo- 
nent a of the power law C{r) oc l/r a decreased from 1.8 
in 1983 to 1.4 in 2000. This means that the upper tail 
became "fatter". Another useful parameter is the total 
income of the upper class as the fraction / of the total 
income in the system. T he fr action / increased from 4% 
in 1983 to 20%' in 2000 [l05{ . However, in year 2001, a 
increased and / decreases, indicating that the upper tail 
was reduced after the stock market crash at that time. 
These results indicate that the upper tail is highly dy- 
namical and not stationary. It tends to swell during the 
stock market boom and shri nk during the bust. Similar 
results were found for Japan fiol [TTa [ilil ITT^ . 

Although relative income inequality within the lower 
class remains stable, the overall income inequality in the 
USA has increased significantly as a result of the tremen- 
dous growth of the income of the upper class. This is 
illustrated by the Lorenz curve and the Gini coefficient 
shown in Fig. [5J The Lorenz curve [99j is a standard way 
of representing income distribution in the economic liter- 
ature. It is defined in terms of two coordinates x(r) and 



y(r) depending on a parameter r: 

r rr'P(r')dr' 

The horizontal coordinate x(r) is the fraction of the pop- 
ulation with income below r, and the vertical coordinate 
y(r) is the fraction of the income this population accounts 
for. As r changes from to oo, x and y change from to 
1 and parametrically defines a curve in the (x, y)-plane. 

Fig. [8J shows the data points for the Loren z curv es in 
1983 and 2000, as computed by the IRS [HI]. Ref - [lQ3| 
analytically derived the Lorenz curve formula y — x + 
(1 — x) ln(l — x) for a purely exponential distribution 
P(r) = cexp(— r/T r ). This formula is shown by the red 
line in Fig.[5]and describes the 1983 data reasonably well. 
However, for year 2000, it is essential to take into account 
the fraction / of income in the up per tail, which modifies 
for the Lorenz formula as follows |l03l . Il04l Il05l | 

y = (1 - f)[x + (1 - x) ln(l -x)]+f Q(x - 1). (20) 

The last term in Eq. (|20[) represent the vertical jump of 
the Lorenz curve at x = 1, where a very small percentage 
of population in the upper class accounts for a substantial 
fraction / of the total income. The blue curve represent- 
ing Eq. (J2DJ) fits the 2000 data in Fig. very well. 

The deviation of the Lorenz curve from the straight 
diagonal line in Fig. [8J is a certain measure of income 
inequality. Indeed, if everybody had the same income, 
the Lorenz curve would be a diagonal line, because the 
fraction of income would be proportional to the fraction 
of the population. The standard measure of income in- 
equality is the so-called Gini coefficient < G < 1, which 
is defined as the area between the Lorenz curve and the 
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FIG. 7 Cumulative probability distribution of tax returns 
plotted on log- log scale versus r/T r (the annual income r nor- 
malized by the average income T r in the exponential part of 
the distribution). The IRS data points are for 1983-2001, and 
the columns of numbers give the values of T r for the corre- 
sponding years. (Reproduced from Ref. [l05| ) 
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FIG. 8 Main panel: Lorenz plots for income distr ibuti on in 
1983 and 2000. The data points are from the IRS [HI, and 
the theoretical curves represent Eq. (|20p with f fr om Fig. 
[7J Inset: The closed circles are the IRS data hjj for the 
Gini coefficient G, and the open circles show t he th eoretical 
formula G = (1 + ,f)/2. (Reproduced from Ref. fw$ ) 
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United States, Bureau of Census data for 1E 
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Annual family income, k$ 

FIG. 9 Histogram: Probability distribution of family income 
for families with two adults (US Census Bureau data). Solid 
line: Fit to Eq. (J2T}. (Reproduced from Ref. [Io3 ].) 



diagonal line, divided by the area of the triangle beneath 
the diagonal line [9^]. Time e volut ion of the Gini coeffi- 
cient, as compu ted b y the IRS [ll3| . is shown in the inset 
of Fig. |H Ref. [loll' derived analytically the result that 
G = 1/2 for a purely exponential distribution. In the 
first approximation, the values of G shown in the inset 
of Fig. [8]are indeed close to the theoretical value 1/2. If 
we take into account the upper tail using Eq. (f2T)]) . the 
form ula for the Gini coefficient becomes G = (1 + /)/2 
|l05j | . The inset in Fig. [8J shows that this formula gives 
a very good fit to the IRS data for the 1990s using the 
values of / deduced from Fig. [71 The values G < 1/2 in 
the 1980s cannot be captured by this formula, because 
the Lorenz data points are slightly above the theoretical 
curve for 1983 in Fig. [51 Overall, we observe that income 
inequality has been increasing for the last 20 years, be- 
cause of swelling of the Pareto tail, but decreased in 2001 
after the stock market crash. 

Thus far we discussed the distribution of individual 
income. An interesting related question is the distribu- 
tion P2{r) of family income r = r\ + where r% and 
r2 are the incomes of spouses. If individual incomes are 
distributed exponentially P(r) oc exp(— r/T r ), then 



P 2 (r) 



dr'P(r')P(r - r') = cr exp(-r/T r ), (21) 



where c is a normalization constant. Fig. [5] shows that 
Eq. (|2ip is in good agreement with the family income 
distribution data from the US Census Bureau [l02j |. In 
Eq. (j2"Tj) , we assumed that incomes of spouses are uncor- 
related. This simple approximation is indeed supported 
by the scatter plot of incomes of spouses shown in Fig. 
1101 Each family is represented in this plot by two points 
{Tx,r2) and (r2,ri) for symmetry. We observe that the 
density of points is approximately constant along the 
lines of constant family income n + ?'2 = const, which 



indicates that incomes of spouses are approximately un- 
corrected. There is no significant clustering of points 
along the diagonal r\ — r 2 , i.e., no strong positive corre- 
lation of spouses' incomes. 

The Gini coefficient for th e fam ily income distribution 
(|2T|) was calculated in Ref. [To2| as G = 3/8 = 37.5%. 
Fig. [11] shows the Lorenz quintiles and the Gini coeffi- 
cient for 1947-1994 plotted from the US Census Bureau 
data. The solid line, representing the Lorenz curve calcu- 
lated from Eq. (j2"Tj) . is in good agreement with the data. 
The systematic deviation for the top 5% of earners results 
from the upper tail, which has a less pronounced effect on 
family income than on individual income, because of in- 
come averaging in the family. The Gini coefficient, shown 
in the inset of Fig. [TTJ is close to the calculated value of 
37.5%. Moreover, the average G for the developed cap- 
italist countries of North America and western Europe, 
as determined by the World Bank 1 1031 ] . is also close to 
the calculated value 37.5%. 

Income distribution has been examined in e cono- 
phy s ics paper s for different countries: J apan l68l. Il09l 
Ilia [ml. Illi ITU JlTsL UTtil , Germany fill El > the 
Ill8j. Italy fllllnallla . the USA 
, Australia [9ll Il20l fl22j . and New 
Zealand [68l [116J. The distributions are qualitatively 
similar to the results presented in this section. The upper 
tail follows a power law and comprises a small fraction of 
population. To fit the lower part of the distribution, dif- 
ferent papers used exponential, Gamma, and log-normal 
distributions. Unfortunately, income distribution is of- 
ten reported by statistical agencies for households, so it 
is difficult to differentiate between one-earner and two- 
earner income distributions. Some papers used interpo- 
lating functions with different asymptotic behavio r fo r 
low and high incomes, such as the Tsallis function lll( 
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FIG. 10 Scatter plot of the spouses' incomes (ri,ra) and 
(r2,ri) based on the data from the Panel Study of Income 
Dynamics (PSID). (Reproduced from Ref. fjQ^ ) 
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United States, Bureau of Census data for 1947-1994 




1 20 30 40 50 60 70 80 90 1 00% 
Cumulative percent of families 

FIG. 11 Main panel: Lorenz plot for family income calculated 
from Eq. (|2ip . compared with the US Census data points. 
Inset: The US Census data points for the Gini coefficient 
for families, compared with the theo retica lly calculated value 
3/8=37.5%. (Reproduced from Ref. (io^ |) 



and the Kaniadakis function [1181 ]. However, the transi- 
tion between the lower and upper classes is not smooth 
for the US data shown in Figs. [6] and [7j so such func- 
tions would not be useful in this case. The special case 
is income distribution in Argentina during the economic 
crisis, whic h sho ws a time-dependent bimodal shape with 
two peaks [l!6j |. 



B. Theoretical models of income distribution 

Having examined the empirical data on income distri- 
bution, let us now discuss theoretical models. Income 
Ti is the influx of money per unit time to an agent i. 
If the money balance rrii is analogous to energy, then 
the income n would be analogous to power, which is the 
energy flux per unit time. So, one should conceptually 
distinguish between the distributions of money and in- 
come. While money is regularly transferred from one 
agent to another in pairwise transactions, it is not typi- 
cal for agents to trade portions of their income. Never- 
theless, indirect transfer of income may occur when one 
employee is promoted and another demoted while the to- 
tal annual budget is fixed, or when one company gets a 
contract whereas another one loses it, etc. A reasonable 
approach, whi c h has a l ong tradition in the economic 
literature jl23l . 1 1241 . Il25l ], is to treat individual income 
r as a stochastic process and study its probability dis- 
tribution. In general, one can study a Markov process 
generated by a matrix of transitions from one income to 
another. In the case where income r changes by a small 
amount Ar over a time period At, the Markov process 
can be treated as income diffusion. Then one can ap- 
ply the general Fokker-Planck equation [7l[ to describe 



evolution in time t of the income distribution function 
P(r,t) EqI 



dP 
~dt 



d_ 

dr 



AP 



d{BP) 



dr 



A 



(Ar) 



B 



{(Ar) 2 ) 



At ' 2At 

(22) 

The coefficients A and B in Eq. (|2"2"|) are determined by 
the first and second moments of income changes per unit 
time. The stationary solution <9 t P = of Eq. (|22|) obeys 
the following equation with the general solution 



d(BP) 
dr 



= -AP, P(r) = 



B(r) 



exp 



Mf) 

B(r') 



dr' 



(23) 

For the lower part of the distribution, it is reasonable 
to assume that Ar is independent of r, i.e., the changes of 
income are independent of i ncom e itself. This process is 
called the additive diffusion [105]. In this case, the coeffi- 
cients in Eq. (|2"2"|) are constants Ag and Bq . Then Eq. (|2"3"|) 
gives the exponential distribution P(r) oc exp(— r/T r ) 
with the effective income temperature T r — Bq/Aq. [No- 
tice that a meaningful stationary solution (|23[) requires 
that A > 0, i.e., (Ar) < 0.] The coincidence of this re- 
sult with the Boltzmann-Gibbs exponential law JTJ) and 
((T]) is not accidental. Indeed, instead of considering pair- 
wise interaction between particles, one can derive Eq. ([T]) 
by considering energy transfers between a particle and a 
big reservoir, as long as the transfer process is "additive" 
and does not involve a Maxwell-demon-like discrimina- 
tion. Stochastic income fluctuations are described by 
a similar process. So, although money and income are 
different concepts, they may have similar distributions, 
because they are governed by similar mathematical prin- 
ciples. It was shown explicitly in Refs. [25|, that 
the models of pairwise money transfer can be described 
in a certain limit by the Fokker-Planck equation. 

On the other hand, for the upper tail of the income 
distribution, it is reasonable to expect that Ar oc r, i.e., 
income changes are proportional to income itself. This 
is known as the proportionality principle of Gibrat 
and the process is called the multiplicative diffusion 
In this case, A = ar and B = br 2 , and Eq. (123[) gives the 
power-law distribution P(r) oc l/r Q+1 with a = 1 + a/b. 

Generally, lower-class income comes from wages and 
salaries, where the additive process is appropriate, 
whereas upper-class income comes from bonuses, in- 
vestments, and capital gains, calculated in p ercentages, 
where the multiplicative process applies [l26| |. However, 
the additive and multiplicative processes may coexist. 
An employee may receive a cost-of-living raise calculated 
in percentages (the multiplicative process) and a merit 
raise calculated in dollars (the additive process). In this 



123], 



105]. 



case, we have A = Ao + ar and B = Bo + br 



b(rl 



where r 2 = Bq jb. Substituting these expressions into Eq. 



we find 



P(r) 



g— (ro/T r ) arctan(r'/r ) 

[1 + {r/r ) 2 ] 1+a / 2b ' 



(24) 
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The distribution (|24[) interpolates between the exponen- 
tial law for low r and the power law for high r, because 
either the additive or the multiplicative process domi- 
nates in the corresponding limit. The crossover between 
the two regimes takes place at r ~ rg, where the addi- 
tive and multiplicative contributions to B are equal. The 
distribution (|24|) has three parameters: the "income tem- 
perature" T r = Aq/Bq, the Pareto exponent a = 1 + a/b, 
and the crossover income r$. It is a minimal model that 
captures the salient features of the empirical income dis- 
tribution shown in Fig. [SJ A mathematically similar, but 
more econ omically oriented model was proposed in Refs. 
[ll4 , Ill5j , where labor income and assets accumulation 
are described by the additive and a multiplicative pro- 
cesses correspondingly. A general stochastic process with 
additive and multiplicative noise was studied numerically 
in Ref. jl27| , but the stationary distribution was not de- 
rived analytically. A similar proc ess w ith discrete time 
increments was studied by Kesten [l28j . Rece ntly, a for- 
mula similar to Eq. (j2"4"|) was obtained in Ref. [129j | . 

To verify the multiplicative and additive hypotheses 
empirically, it is necessary to have data on income mo- 
bility, i.e., the income changes Ar of the same people 
from one year to another. The distribution of income 
changes P(Ar\r) conditional on income r is generally not 
available publicly, although it can be reconstructed by 
researchers at the tax agencies. Nevertheless, the mul- 
tiplicative hypothesis f or th e upp er class was quantita- 
tively verified in Refs. [llll Ill2j | for Japan, where tax 
identification data are published for the top taxpayers. 

The power-law distribution is meaningful only when 
it is limited to high enough incomes r > tq. If all 
incomes r from to oo follow a purely multiplicative 
process, then one can change to a logarithmic variable 
x = ln(r/r») in Eq. (|2"2"|) and show that it gives a Gaus- 
sian time-dependent distribution Pt(x) oc cxp(— x 2 /2<r 2 t) 
for x, i.e., the log- normal distribution for r, also known as 
the Gibrat distribution [123]. However, the width of this 
distribution increases linearly in time, so the distrib ution 
is not stationary. This was pointed out by Kalecki [l24j 
a long time ago, but the log-normal distribution is still 
widely used for fitting income distribution, despite this 
fundamental logical flaw in its justification. In the classic 
paper [125 ] , Champernowne showed that a multiplicative 
process gives a stationary power-law distribution when a 
boundary condition is imposed at r$ ^ . Later, this 
result was rediscovered by econophysicists jl30l . Il3l| . In 
our Eq. the exponential distribution of the lower 

class effectively provides such a boundary condition for 
the power law of the upper class. Notice also that Eq. 
([24]) reduces to Eq. (| 18[) in the limit ro — > 0, which cor- 
responds to purely multiplicative noise B = br 2 . 

There are alternative approaches to income distribu- 
tion in econ omic literature. One of them, proposed by 
Lydall jl32| . involves social hierarchy. Groups of peo- 
ple have leaders, which have leaders of a higher order, 
and so on. The number of people decreases geometri- 
cally (exponentially) with the increase of the hierarchical 



level. If individual income increases by a certain factor 
(i.e., multiplicatively) when moving to the next hierar- 
chical level, then income distribution follows a power law 
|l32j . However, the original argument of Lydall can be 
easily modified to produce the exponential distribution. 
If individual income increases by a certain amount, i.e., 
income increases linearly with the hierarchical level, then 
income distribution is exponential. The latter process 
seems to be more realistic for moderate annual incomes 
below 100 k$. A similar scenario is the Bernoulli trials 
[l33j . where individuals have a constant probability of 
increasing their income by a fixed amount. We see that 
the deterministic hierarchical models and the stochas- 
tic models of additive and multiplicative income mobility 
represent essentially the same ideas. 



VI. OTHER APPLICATIONS OF STATISTICAL PHYSICS 

Statistical physics was applied to a number of other 
subjects in economics. Because of lack of space, only two 
such topics are briefly discussed in this section. 



A. Economic temperatures in different countries 

As discussed in Sees. IIV.CI and IV. Al the distributions 
of money, wealth, and income are often described by ex- 
ponential functions for the majority of the population. 
These exponential distributions are characterized by the 
parameters T m , T w , and T r , which are mathematically 
analogous to temperature in the Boltzmann-Gibbs dis- 
tribution |T]). The values of these parameters, extracted 
from the fits of the empirical data, are generally different 
for different countries, i.e., different countries have dif- 
ferent economic "temperatures". For example, Ref. (95j 
found that the US income temperature was 1.9 times 
higher than the UK income temperature in 1998 (using 
the exchange rate of dollars to pounds at that time). 
Also, there was ±25% variation between income temper- 
atures of different states within the USA [95| . 

In physics, a difference of temperatures allows one to 
set up a thermal machine. In was argued in Ref. [25[ that 
the difference of money or income temperatures between 
different countries allows one to extract profit in inter- 
national trade. Indeed, as discussed at the end of Sec. 
IIV.CI the prices of goods should be commensurate with 
money or income temperature, because otherwise people 
cannot afford to buy those goods. So, an international 
trading company can buy goods at a low price T\ in a 
"low-temperature" country and sell them at a high price 
T2 in a "high-temperature" country. The difference of 
prices T2 — T1 would be the profit of the trading company. 
In this process, money (the analog of energy) flows from 
the "high-temperature" to the "low-temperature" coun- 
try, in agreement with the second law of thermodynam- 
ics, whereas products flow in the opposite direction. This 
process very much resembles what is going on in global 



17 



economy now. In this framework, the perpetual trade 
deficit of the USA is the consequence of the second law 
of thermodynamics and the difference of temperatures 
between the USA and the "low-temperature" countries, 
such as Chin a. Similar ideas were developed in more de- 
tail in Refs. [l34l Il35l | , including a formal Carnot cycle 
for international trade. 

The statistical physics approach demonstrates that 
profit originates from statistical nonequilibrium (the dif- 
ference of temperatures), which exists in the global econ- 
omy. However, it does not answer the question what is 
the origin of this difference. By analogy with physics, one 
would expect that the money flow should reduce the tem- 
perature difference and, eventually, lead to equilibriza- 
tion of temperatures. In physics, this situation is known 
as the "thermal death of the universe" . In a completely 
equilibrated global economy, it would be impossible to 
make profit by exploiting differences of economic tem- 
peratures between different countries. Although global- 
ization of modern economy does show a tendency toward 
cquilibrization of living standards in different countries, 
this process is far from straightforward, and there are 
many examples contrary to equilibrization. This interest- 
ing and timely subject certainly requires further study. 



cists in dealing with phase transitions and alloys may be 
helpful for practical applications of such models [l42j | . 

VII. FUTURE DIRECTIONS, CRITICISM, AND 
CONCLUSIONS 

The statistical models described in this review are 
quite simple. It is commonly accepted in physics that 
theoretical models are not intended to be photographic 
copies of reality, but rather be caricatures, capturing the 
most essential features of a phenomenon with a minimal 
number of details. With only few rules and parameters, 
the models discussed in Sees. IIIIl IIV1 and [V] reproduce 
spontaneous development of stable inequality, which is 
present in virtually all societies. It is amazing that the 
calculated Gini coefficients, G = 1/2 for individuals and 
G = 3/8 for families, are actually very close to the US 
income data, as shown in Fig. [5] and QT] These simple 
models establish a baseline and a reference point for de- 
velopment of more sophisticated and more realistic mod- 
els. Some of these future directions are outlined below. 



A. Future directions 



B. Society as a binary alloy 

In 1971, Thomas Schelling propo sed the now-famous 
mathematical model of segregation jl36l |. He considered 
a lattice, where the sites can be occupied by agents of two 
types, e.g., blacks and whites in the problem of racial seg- 
regation. He showed that, if the agents have some prob- 
abilistic preference for the neighbors of the same type, 
the system spontaneously segregates into black and white 
neighborhoods. This mathematical model is similar to 
the so-called Ising model, which is a popular model for 
studying phase transitions in physics. In this model, each 
lattice site is occupied by a magnetic atom, whose mag- 
netic moment has only two possible orientations, up or 
down. The interaction energy between two neighboring 
atoms depends on whether their magnetic moments point 
in the same or in the opposite directions. In physics lan- 
guage, the segregation found by Schelling represents a 
phase transition in this system. 

Another similar model is the binary alloy, a mixture of 
two elements w hich attract or repel each other. It was 
noticed in Ref. [l37l | that the behavior of actual binary 
alloys is striking ly simila r to social segregation. In the 
following papers [42|, Il38| , this mathematical analogy was 
developed further and compared with social data. Inter- 
esting concepts, such as the coexistence curve between 
two phases and the solubility limit, were discussed in this 
work. The latter concept means that a small amount of 
one substance dissolves into another up to some limit, 
but phase separation (segregation) develops for higher 
concentr ations. Recentl y, similar ideas were rediscovered 
m Refs. [i"39Lfl40l. Il4l| . The vast experience of physi- 



a. Agents with a finite lifespan. The models discussed in 
this review consider immortal agents who live forever, like 
atoms. However, humans have a finite lifespan. They 
enter the economy as young people and exit at an old 
age. Evolution of income and wealth as functions of age 
is studied in economics using the so-called overlapping- 
generations model. The absence of the age variable was 
one of the criticisms of econophysics by the economist 
Paul Anglin [3l[. However, the drawback of the stan- 
dard overlapping-generations model is that there is no 
variation of income and wealth between agents of the 
same age, because it is a representative-agent model. It 
would be best to combine stochastic models with the age 
variable. Also, to take into account inflation of average 
income, Eq. (|22[) should be rewritten for relative income, 
in the spirit of Eq. (fl~7|) . These modifications would allow 
to study the effects of demographic waves, such as baby 
boomers, on the distributions of income and wealth. 

b. Agent-based simulations of the two-class society. The 
empirical data presented in Sec. lV.Al show quite convinc- 
ingly that the US population consists of two very distinct 
classes characterized by different distribution functions. 
However, the theoretical models discussed in Secs. lIIII and 
HVI do not produce two classes, although they do produce 
broad distributions. Generally, not much attention has 
been payed in the agent-based literature to simul ation 
of two classes. One important exception is Ref. [l43j |. 
in which spontaneous development of employers and em- 
ployees classes from initially equal agents was simulated 
[36| , More work in this direction would be certainly de- 
sirable. 
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c. Access to detailed empirical data. A great amount of 
statistical information is publicly available on the Inter- 
net, but not for all types of data. As discussed in Sec. 
IIV.C1 it would be very interesting to obtain data on the 
distribution of balances on bank accounts, which would 
give information about the distribution of money (as op- 
posed to wealth). As discussed in Sec. IV. Bl it would be 
useful to obtain detailed data on income mobility, to ver- 
ify the additive and multiplicative hypotheses for income 
dynamics. Income distribution is often reported as a mix 
of data on individual income and family income, when 
the counting unit is a tax return (joint or single) or a 
household. To have a meaningful comparison with theo- 
retical models, it is desirable to obtain clean data where 
the counting unit is an individual. Direct collaboration 
with statistical agencies would be very useful. 



d. Economies in transition. Inequality in developed capi- 
talist countries is generally quite stable. The situation is 
very different for the former socialist countries making a 
transition t o a m arket economy. According to the World 
Bank data jl03l |. the average Gini coefficient for family 
income in eastern Europe and the former Soviet Union 
jumped from 25% in 1988 to 47% in 1993. The Gini 
coefficient in the socialist countries before the transition 
was well below the equilibrium value of 37.5% for mar- 
ket economies. However, the fast collapse of socialism 
left these countries out of market equilibrium and gen- 
erated a much higher inequality. One may expect that, 
with time, their inequality will decrease to the equilib- 
rium value of 37.5%. It would be very interesting to trace 
how fast this relaxation takes place. Such a study would 
also verify whether the equilibrium level of inequality is 
universal for all market economies. 



e. Relation to physical energy. The analogy between en- 
ergy and money discussed in Sec. IIII.Bl is a formal math- 
ematical analogy. However, actual physical energy with 
low entropy (typically in the form of fossil fuel) also plays 
a very important role in the modern economy, being the 
basis of current human technology. In view of the loom- 
ing energy and climate crisis, it is imperative to find real- 
istic ways for making a transition from the current "dis- 
posable" economy based on "cheap" and "unlimited" en- 
ergy and natural resources to a sustainable one. Hetero- 
geneity of human society is one of the important factors 
affecting such a transition. Econophysics, at the intersec- 
tion of energy, entropy, economy, an d sta tistical physics, 
may play a useful role in this quest |l44j . 



B. Criticism from economists 

As econophysics is gaining popularity, some criticism 
has appeared from economists [3l|, including those who 
are closely involved with the econophysics movement 



[13, [33, 34]. This reflects a long-standing tradition in eco- 
nomic and social sciences of writing critiques on different 
schools of thought. Much of the criticism is useful and 
constructive and is already being accommodated in the 
econophysics work. However, some criticism results from 
misunderstanding or miscommunication between the two 
fields and some from significant differences in scientific 
philosophy. Several insightful respo nses to the critic ism 
have been published |l45l . Il46l [l47j . see also [7l.ll48|. In 
this section, we briefly address the issues that are directly 
related to the material discussed in this review. 



a. Awareness of pre vio us e cono mic literature. One com- 
plaint of Refs. [3l|, [H, HH is that physicists are not 
well aware of the previous economic literature and ei- 
ther rediscover known results or ignore well-established 
approaches. To address this issue, it is useful to keep 
in mind that science itself is a complex system, and sci- 
entific progress is an evolutionary process with natural 
selection. The sea of scientific literature is enormous, 
and nobody knows it all. Recurrent rediscovery of reg- 
ularities in the natural and social world only confirms 
their validity. Independent rediscovery usually brings a 
different perspective, broader applicability range, higher 
accuracy, and better mathematical treatment, so there is 
progress even when some overlap with previous results 
exists. Physicists are grateful to economists for bringing 
relevant and specific references to their attention. Since 
the beginning of modern econophysics, many old refer- 
ences have been uncovered and are now routinely cited. 

However, not all old references are relevant to the new 
development. For example, Ref. [33| complained that the 
econophysics literature on income distribution ignores 
the so-called Kuznets hypothesis [l49l | . The Kuznets hy- 
pothesis postulates that income inequality first rises dur- 
ing an industrial revolution and then decreases, produc- 
ing an inverted-U-shaped curve. Ref. [H| admitted that, 
to date, the large amount of literature on the Kuznets 
hypothesis is inconclusive. Ref. [HJ mentioned that this 
hypothesis applies to the period from colonial times to 
1970s; however, the empirical data for this period are 
sparse and not very reliable. The econophysics literature 
deals with the reliable volumes of data for the second half 
of the 20th century, collected with the introduction of 
computers. It is not clear what is the definition of indus- 
trial revolution and when exactly it starts and ends. The 
chain of technological progress seems to be continuous 
(steam engine, internal combustion engine, cars, plastics, 
computers, Internet), so it is not clear where the pur- 
ported U-curve is supposed to be placed in time. Thus, 
the Kuznets hypothesis appears to be, in principle, unver- 
ifiable and unfalsifiable. The original paper by Kuznets 
[l49j | actually does not contain any curves, but it has 
one table filled with made-up, imaginary data! Kuznets 
admits that he has "neither the neces sary data nor a rea- 
sonably complete theoretical model" 149, p 12]. So, this 
paper is understandably ignored by the econophysics lit- 
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erature. In fact, the data analysis for 1 947- 1984 shows 
amazing stability of income distribution |l50l ] , consistent 
with Fig. [TlJ The increase of inequality in the 1990s re- 
sulted from growth of the upper tail relative to the lower 
class, but the relative inequality within the lower class 
remains very stable, as shown in Fig. [7] 



b. Reliance on visual data analysis. Another complaint, of 
Ref. [33l | is that econophysicists favor graphic analysis of 
data over the formal and "rigorous" testing prescribed by 
mathematical statistics, as favored by economists. This 
complaint goes against the trend of all sciences to use in- 
creasingly sophisticated data visualization for uncovering 
regularities in complex system. The thick IRS publica- 
tion 1304 [15 it ] is filled with data tables, but has virtu- 
ally no graphs. Despite the abundance of data, it gives a 
reader no idea about income distribution, whereas plot- 
ting the data immediately gives insight. However, in- 
telligent plotting is the art with many tools, which not 
many researchers have mastered. The author completely 
agrees with Ref. [33| that too many papers mindlessly 
plot any kind of data on a log-log scale, pick a finite 
interval, where any smooth curved line can be approxi- 
mated by a straight line, and claim that there is a power 
law. In many cases, replotting the same data on a log- 
linear scale converts a curved line into a straight line, 
which means that the law is actually exponential. 

Good visualization is extremely helpful in identifying 
trends in complex data, which can then be fitted to a 
mathematical function. However, for a complex system, 
such a fit should not be expected with infinite precision. 
The fundamental laws of physics, such as Newton's law 
of gravity or Maxwell's equations, are valid with enor- 
mous precision. However, the laws in condensed mat- 
ter physics, uncovered by experimentalists with a combi- 
nation of visual analysis and fitting, usually have much 
lower precision, at best 10% or so. Most of these laws 
would fail the formal criteria of mathematical statistics. 
Nevertheless these approximate laws are enormously use- 
ful in practice, and the everyday devices, engineered on 
the basis of these laws, work very well for all of us. 

Because of the finite accuracy, different functions may 
produce equally good fits. Discrimination between the 
exponential, Gam ma, a nd log-normal functions may not 
be always possible [122]. However, the exponential func- 
tion has fewer fitting parameters, so it is preferable on 
the basis of simplicity. The other two functions can 
simply mimic the exponential function with a particu- 
lar choice of the additional parameters [l22j. Unfortu- 
nately, many papers in mathematical statistics introduce 
too many fitting parameters into complicated functions, 
such as the generalized beta distribution mentioned in 
Ref. [33] ]. Such overparametrization is more misleading 
than insightful for data fitting. 



c. Quest for universality. Ref. [33| criticized physicists for 
trying to find universality in economic data. It also seems 
to equate the concepts of power law, scaling, and uni- 
versality. These are three different, albeit overlapping, 
concepts. Power laws usually apply only to a small frac- 
tion of data at the high ends of various distributions. 
Moreover, the exponents of these power laws are usually 
nonuniversal and vary from case to case. Scaling means 
that the shape of a function remains the same when its 
scale changes. However, the scaling function does not 
have to be a power-law function. A good example of 
scaling is shown in Fig. [7] where income distributions for 
the lower class collapse on the same exponential line for 
about 20 years of data. We observe amazing universal- 
ity of income distribution, unrelated to a power law. In a 
general sense, the diffusion equation is universal, because 
it describes a wide range of systems, from dissolution of 
sugar in water to a random walk in the stock market. 

Universalities are not easy to uncover, but they form 
the backbone of regularities in the world around us. This 
is why physicists are so much interested in them. Uni- 
versalities establish the first-order effect, and deviations 
represent the second-order effect. Different countries may 
have somewhat different distributions, and economists of- 
ten tend to focus on these differences. However, this focus 
on details misses the big picture that, in the first approx- 
imation, the distributions are quite similar and universal. 



d. Theoretical modeling of money, wealth, and income. 

Refs. [3ll. I33l.l34| pointed out that many econophysics pa- 
pers confuse or misuse the terms for money, wealth, and 
income. It is true that terminology is sloppy in many pa- 
pers and should be refined. However, the terms in Refs. 
[25l l2ll | are quite precise, and this review clearly distin- 
guishes between these concepts in Sees. Mil IIV1 andfVl 

One contentious issue is about conservation of money. 
Ref. [33| agrees that "transactions are a key economic 
process, and they are necessarily conservative", i.e., 
money is indeed conserved in transactions between 
agents. However, Refs. jU [34| complain that the 
models of conservative exchange do not consider produc- 
tion of goods, which is the core economic process and 
the source of economic growth. Material production is 
indeed the ultimate goal of the economy, but it does not 
violate conservation of money by itself. One can grow cof- 
fee beans, but nobody can grow money on a money tree. 
Money is an artificial economic device that is designed to 
be conserved. As explained in Sec. IIIIl the money trans- 
fer models implicitly assume that money in transactions 
is voluntarily payed for goods and services generated by 
production for the mutual benefit of the parties. In prin- 
ciple, one can introduce a billion of variables to keep track 
of every coffee bean and other product of the economy. 
What difference would it make for the distribution of 
money! Despite claims in Refs. (3lL l33l]. there is no con- 
tradiction between models of conservative exchange and 
the classic work of Adam Smith and David Ricardo. The 
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difference is oniy in the focus: We keep track of money, 
whereas they keep track of coffee beans, from produc- 
tion to consumption. These approaches address different 
questions, but do not contradict each other. Because 
money constantly circulates in the system as payment 
for production and consumption, the resulting statistical 
distribution of money may very well not depend on what 
is exactly produced and in what quantities. 

In principle, the models with random transfers of 
money should be considered as a reference point for de- 
veloping more sophisticated models. Despite the totally 
random rules and "zero intelligence" of the agents, these 
models develop well-characterized, stable and stationary 
distributions of money. One can modify the rules to 
make the agents more intelligent and realistic and see 
how much the resulting distribution changes relative to 
the reference one. Such an attempt was made in Ref. 
[32I ] by modifying the model of Ref. [77| with various 
more realistic economic ingredients. However, despite 
the modifications, the resulting distributions were essen- 
tially the same as in the original model. This example 
illustrates the typical robustness and universality of sta- 
tistical models: Modifying details of microscopic rules 
does not necessarily change the statistical outcome. 

Another misconception, elaborated in Ref. [13, S| , is 
that the money transfer models discussed in Sec. Illll im- 
ply that money is transferred by fraud, theft, and vio- 
lence, rather than voluntarily. One should keep in mind 
that the catchy labels "theft-and-fraud" , "marriage-and- 
divorce" , and "yard-sale" were given to the money trans- 
fer mo dels by the journalist Brian Hayes in a popular 
article [l52j |. Econophysicists who originally introduced 
and studied these models do not subscribe to this termi- 
nology, although the early work of Angle [5l[ did mention 
violence as one source of redistribution. In the opinion 
of the author, it is indeed difficult to justify the propor- 
tionality rule ([5]), which implies that agents with high 
balances pay proportionally greater amounts in transac- 
tions than agents with low balances. However, the addi- 
tive model of Ref. [25| . where money transfers Am are 
independent of money balances m, of the agents, docs 
not have this problem. As explained in Sec. IIII.C1 this 
model simply means that all agents pay the same prices 
for the same product, although prices may be different 
for different products. So, this model is consistent with 
voluntary transactions in a free market. 

Ref. [l45( | argued that conservation of money is vio- 
lated by credit. As explained in Sec. IIII.D1 credit does 
not violate conservation law, but creates positive and 
negative money without changing net worth. Neg ative 
money (debt) is as real as positive money. Ref. [l45| | 
claimed that money can be easily created with the tap of 
a computer key via credit. Then why would an employer 
not tap the key and double salaries, or a funding agency 
double research grants? Because budget constraints are 
real. Credit may provide a temporary relief, but sooner 
or later it has to be paid back. Allowing debt may pro- 
duce a double-exponential distribution as shown in Fig. 



|3l but it does not change the distribution fundamentally. 

As discussed in Sec. IIII.B1 a central bank or a central 
government can inject new money into the economy. As 
discussed in Sec. IIV1 wealth is generally not conserved. 
As discussed in Sec. [Vj income is different from money 
and is described by a different model ([22]). However, 
the empirical distribution of income shown in Fig. [5] is 
qualitatively similar to the distribution of wealth shown 
in Fig. and we do not have data on money distribution. 

C. Conclusions 

The "invasion" of physicists into economics and fi- 
nance at the turn of the millennium is a fascinating 
phenomenon. The physicist Joseph McCauley proclaims 
that "Econophysics will displace economics in both the 
universities and boardrooms, simply becau se w hat is 
taught in economics classes doesn't work" Il53l. Al- 
though there is some truth in his arguments [1451 ] . one 
may consider a less radical scenario. Econophysics may 
become a branch of economics, in the same way as games 
theory, psychological economics, and now agent-based 
modeling became branches of economics. These branches 
have their own interests, methods, philosophy, and jour- 
nals. The main contribution from the infusion of new 
ideas from a different field is not in answering old ques- 
tions, but in raising new questions. Much of the misun- 
derstanding between economists and physicists happens 
not because they are getting different answers, but be- 
cause they are answering different questions. 

The subject of income and wealth distributions and 
social inequality was very popular at the turn of an- 
other century and is associated with the names of Pareto, 
Lorenz, Gini, Gibrat, and Champernowne, among others. 
Following the work by Pareto, attention of researchers 
was primarily focused on the power laws. However, when 
physicists took a fresh, unbiased look at the empirical 
data, they found a different, exponential law for the lower 
part of the distribution. The motivation for looking at 
the exponential law, of course, came from the Boltzmann- 
Gibbs distribution in physics. Further studies provided 
a more detailed picture of the two-class distribution in a 
society. Although social classes have been known in po- 
litical economy since Karl Marx, realization that they are 
described by simple mathematical distributions is quite 
new. Demonstration of the ubiquitous nature of the ex- 
ponential distribution for money, wealth, and income is 
one of the new contributions produced by econophysics. 
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